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Abstract 

We study here a heat-type differential equation of order n greater than 
two, in the case where the time-derivative is supposed to be fractional. 
The corresponding solution can be described as the transition function of 
a pseudoprocess ^„ (coinciding with the one governed by the standard, 
non-fractional, equation) with a time argument Ta which is itself random. 
The distribution of Ta is presented together with some features of the 
solution (such as analytic expressions for its moments). 
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1 Introduction 

The study of diffusion equations with a fractional derivative component have 
been firstly motivated by the analysis of thermal diffusion in fractal media in 
NigmatuUin (1986) and Saichev and Zaslavsky (1997). This topic has been 
extensively treated in the probabilistic literature since the end of the Eighties: 
see, for examples, Wyss (1986), Schneider and Wyss (1989), Mainardi (1996), 
Angulo et al. (2000). Recently fractional equations of different types have been 
also studied, such as, for example, the Black and Scholcs equation (see Wyss 
(2000)) and the fractional diffusion equations with stochastic initial conditions 
(see Anh and Leonenko (2000)). 

Our aim will concern the extension, to the case of fractional time-derivative, 
of a class of equations which is well known in the literature, namely the higher- 
order heat-type equations. Therefore we will be interested in the solution of the 
following problem, for < a < 1, n > 2, 

^u{x,t) ^ kn-§^u[x,t) X e M, < > 0, ...s 
u{x,0) = 5{x) 
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where 5{-) is the Dirac delta function, fc„ = (—1)'"'"-'^ for n — 2q, q £ N, while 
fc„ = ±1 for n — 2q + 1. The fractional derivative appearing in ([1]) is meant, in 
the Dzherbashyan-Caputo sense, as 



where m — I — [aj and / e C™ (see Samko et al. (1993) for a general reference 
on fractional calculus). 

In the non-fractional case (which can be obtained from ours as a particular 
case, for a = 1) the pseudoprocesses ^'^ — 5'n(i), t > driven by n-th order 
equations, i.e. 



for n > 2, have been introduced in the Sixties and studied so far by many authors 
starting from Krylov (1960), Daletsky (1969). Moreover the distributions of 
many functionals of \E'„ have been obtained: in Hochberg and Orsingher (1994) 
the distribution of sojourn time on the positive half-line is presented, for n odd, 
while for an arbitrary n the same topic is analyzed in Lachal (2003). For n = 3, 4, 
the case where the pseudoprocess is constrained to be zero at the end of the time 
interval is considered in Nikitin and Orsingher (2000) and the corresponding 
distribution of the sojourn time is evaluated. In Beghin et al. (2000) the 
distribution of the maximum is obtained under the same circumstances. In the 
unconditional case the maximal distribution is presented in Orsingher (1991), 
for n odd, while the joint distribution of the maximum and the process for 
diffusion of order n = 3,4 is presented in Beghin et al. (2001). Lachal (2003) 
has extended these results to any order n > 2. 

Some other functionals, such as the first passage time, are treated in Nishioka 
(1997) and Lachal (2007). Finally in Beghin and Orsingher (2005) it is proved 
that the local time in zero possesses a proper probability distribution which 
coincides with the (folded) solution of a fractional diffusion equation of order 
2{n- l)/n, n>2. 

In the fractional case under investigation (i.e. for < a < 1) we prove that 
the process related to ([T]) is a pseudoprocess evaluated at a random time 
Ta — Ta{t),t > 0, so that we can write it as 5'„(X^). The probability law of the 
random time is shown to solve a fractional diffusion equation of order 2a and 
it can be expressed in terms of Wright functions. It is interesting to stress that 
the introduction of a fractional time-derivative exerts its influence only on the 
"temporal" argument, while the governing process is not affected and depends 
only on the degree n of the equation. 

Moreover, in section 2, some particular cases of these results are analyzed: 
in the non-fractional case, a = 1, we easily get Ta{t) t. For a = 1/2 it can 
be verified that 7^ coincides with the reflecting Brownian motion and then the 
pseudoprocess governed by equation ([T]) reduces to 'I'„(|i?(i)|), t > (where B 
denotes a standard Brownian motion). 




for m — 1 < a < m 




■p{x,t), 



X e M, t > 0, 



(2) 
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In section 3 the moments of ^n{Ta) are obtained in two alternative ways. 

Section 4 presents some more explicit forms of the solution to equation ([1]) 
which can be obtained by splitting the interval of values for a in two different 
ones and treating them separately. Indeed, if we restrict ourselves to the case 
a G [1/2, 1], the distribution of the random time Ta coincides with ^p2_{u;t), 
u > 0, where by pi. (•; t) we have denoted a stable law of index 1/a. 

As far as the other interval is concerned (i.e. a E (0, 1/2]), an explicit expres- 
sion of the solution can be evaluated by specifying a — 1/m, m S N, m > 2. In 
this particular setting the pseudoprocess can be represented by (G(t)) , t > 0, 
where G{t) — H^Ti Gj{t),t > and Gj{t), j = l,...,m— 1 are independent 
random variables whose density is presented in explicit form. 

Finally we obtain some interesting results by specifying ([1]) for particular 
values of n. For example, taking n = 2, we can conclude that the process related 
to the fractional diffusion equation 

^uix, t) = ^u{x, t) X e M, i > 0, (3) 

for < a < 1, is represented by B(Ta), in accordance with the results already 
known on ([3]). In particular for a — 1/m, equation ^ turns out to be solved 
by the density of the process B {G{t)) ,t > 0. 

In the special case n = 3 the results above reduces to those presented in De 
Gregorio (2002), while, for n = 4, they represent a probabilistic alternative to 
the analytic approach provided by Agrawal (2000). 



2 First expressions for the solution 

We start by considering the n-th order fractional equation and the following 
corresponding initial- value problem, for 0<a<l, n>2, 

w(x,0)=(5(x) ^' 

where fc„ = (— 1)"^+^ for n = 2g, g S N, while kn = ±1 for n = 2(7 -I- 1 and 5{-) 
is the Dirac delta function. The first step consists in evaluating the Laplace 
transform of the solution Ua{x, t), namely 

/"OO 

Uaix,s)= / e^''*Uaix,t)dt, (5) 

and recognizing that it is related to the Laplace transform of the solution p„ (x, t) 
of the corresponding non-fractional n-th order equation ([2]) (which can be de- 
rived from 111 for a = 1). The latter is usually expressed as 

1 /■+°° 

Pu{x,t)^— e^--+'="*(^-)"dz. (6) 

27r . 1^, 
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Theorem 2.1 Let $„(a;, s) — /g°° e ^*-pn{x^t)dt be the Laplace transform of 
the solution to ([2]); then (O can be expressed as follows 



Ua(x,s) 



(7) 



Proof 

By taking the Laplace transform of ^ and considering the initial condition, 
we get 

s^U^{x, s) - = k^^U^{x, s). (8) 

ox" 

Then, by integrating ([5]) with respect to x in [—£,£] and letting e ^ 0, we 
have the following condition for the (n — l)-th derivative 



dx 



— Ua{x,s) 



which must be added to the continuity conditions in zero holding for the j-th 
derivatives, for j = 0,..,n — 2. Therefore our problem is reduced to the n-th 
order linear equations 



fc„^J7a(a;,s) = s"f7a(x, s), x^O 

x=0+ 



dxi ^a{x, S 



x=Q- 
x=0+ 



x=0- 



0, for j = 0, l,...,n-2 



(9) 



If we now impose the boundedness condition for x — )■ ±oo, we obtain 



Ua{x,s) 



if X > 
if X < 



(10) 



where 6k are the n-th roots of fc„, / = {fc : Re(0fe) < 0} and J = {k : Me(6'fe) > 0} . 
The n unknown constants Cfe,fc e / and dk,k £ J, appearing in (jlOp must be 
determined by taking into account the matching conditions in ([9]), as follows: 



Efce/ cfc^i - Efee,7 dkOi = 0, for J = 0, n - 2 
By defining 



(11) 



Cfc, 

-dk, 



like I 
like J ' 



Zk = ■{ J :-f L ^ 7 > (12) 

the linear system in (|lip can be rewritten as the following Vandermonde system 



A:=0 



0, for j = 0,...,n-2 

-fc„s"/"-i, forj = n-l 



(13) 
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Following an argument similar to Beghin and Orsinghcr (2005) (see p. 1024-5) 
we get 



Zk 



n— 1 _j 

= [-irkr^S-h^-'W—— (14) 



r=0 



j -iW"-ie^, iffc„ = l 

\ -iW"-ie^-^, iffc„ = -l 

where, in the last step, we have used formula (2.19) obtained therein. We now 
substitute into (fTU]) the constants evaluated in p^ . taking into account p2p 
and distinguishing the case of n even from the odd one. Indeed, for n = 2(7 + 1, 
the roots of kn are respectively 

J e^, for /c„ = 1 
0.= IIMI)^ (15) 

It- , iUi ft/^ — -L 

SO that (fTO|) becomes, in this case, 

= 1 7 "T'T' Y:""' ^ > « . (16) 

\ lW"-iE,g,,f?fee«'=^ ' fora;<0 

, , I -, (2fc + g + l)7rz 

Analogously, for n~2q and fc„ = (— 1)'+ , the roots are 9k — e " so 
that we get 

( {2k + q+l)TTi 2fcxi 



e " = e 1 , for fc„ = 1 

(2fc + <;+l),ri (2fc+l),ri 

e " = e " , tor k„ = — 1 



(17) 



where, in the first line, we have used the following relationship 

e(«+i)" = (-1)'+^ = fc„ = 1, 
while, in the second one, we have considered the fact that 

= (-l)fc„ = 1. 

Since (fT7)) coincides with we obtain even for n — 2q formula (fTC|) . The 
proof is completed by comparing it with formula (12) of Lachal (2003), which 
reads 

-^s'^""'Efc6/^fce'''=^''"^ for:E>0 
^s'/""'Efeej^fee^''''""> forx<0 

□ 



$„(x,s) 



By inverting the Laplace transform ([7]) we can obtain a first expression of 
the solution in terms of a fractional integral of a particular stable law. Following 
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the notation of Samorodnitsky and Taqqu (1994), we will denote by Sa{<J, /?, A*) 
the distribution of a stable random variable X of index a, with characteristic 
function 

Ee'"^ =exp^^-a"'\s\" (l-iP{sign s) tan— ^ +ifisY a^l, seR. 

^ (18) 
Moreover let I(i~a) denote the Riemann-Liouville fractional integral of order 
1 — a, which is defined as I{i-a) [/(w)] {t) = p^^^^^ J^{t — w)^" f{w)dw (see 
Samko et al. (1993), p.33). 

Theorem 2.2 Let p^{-;u) be the stable distribution 5*0, (cr, 1,0), with param- 
eters a = (ucos7ra/2)^/",/3 — 1, fi = 0, then the fundamental solution to ([4]) 
can be expressed, for < a < 1, as 



Ua{x,t) = p„{x,u)I(_i_a)[Pa{w;u)]{t)du. (19) 
Jo 

Proof 

We recall that, for < a < 2 and a 1, a. stable random variable X ~ 
Sa{(J, 1,0) has Laplace transform 

Eie-"^) = e"-=('°/2)''° , s > 

(see Samorodnitsky and Taqqu (1994), p. 15, for details), so that, in our case 
(for cr = (MCOS7ra/2)i/"), it reduces to E^e^''^) 
rewrite (O as 



OO / /' + 00 



C/„(x,s) = s°'-' / e~" 'pn{x,t)dt (20) 

' _ _ \ 

e ''^^p^{z;u)dz \ pn{x,u)du 
s""^ / e"*^ ( / Pa{z;u)pn{x,u)du\ dz. 



f+OO / /' + 00 



For < a < 1 the first term appearing in ((20l) can be easily inverted by 
considering that 



1 



3°'-^ = — / e-'H-'^dt 

r(i - a) io 

SO that the inverse Laplace transform of (PD)) can be written as 



+ 00 



1 



t / /'+00 



Ua{x,t) = r(l - a) J "( ' Paiw,u)pnix,u)du] dw (21) 



1 



r(i - a) Jo 



+00 

(f — w)^°'p^{w; u)dw ) p„(a;, u)du. 
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Finally we recognize in the last expression a fractional Riemann-Liouville 
integral I{i-a) of order 1 — a of the stable density (where the integration is in- 
tended with respect to the first argument, since the second represents a constant 
in the scale parameter). □ 

The previous result suggests that the solution to our problem can be de- 
scribed as the transition function p„ = pn(x,u) of a pseudoprocess Vl/^ with 
a time-argument 7^ which is itself random. Only for a = 1 we can derive 
from Theorem 2.1 the obvious result that Ta{t) i, so that the solution to 
^ coincides, as expected, with pn{x,t). In all other cases the governing pro- 
cess coincides with the non-fractional one, while the introduction of a fractional 
time-derivative exerts its influence only on the time argument (as remarked 
above) . 

To check that 7^ possesses a true probability density we can observe that 
it is non-negative: indeed it coincides with the fractional integral of a stable 
density S'q((t, 1,0) with skewness parameter equal to 1 (which, by the way, for 
< a < 1, has support restricted to [0, cxd)). Moreover it integrates to one, as 
can be ascertained by the following steps: 



du 



r(i-a)yo 



(t — w) '^p^{w;u)dw 



Til -a) 


Jo 


1 




27rir(l - 


a) J 


1 




27rir(l - 


a) J 


1 




r(a)r(i - 


-a) 



-\-ioo 



' du / {t-wydw — : / e'^^e-' ""ds 
{t - wydw 

] Jo J -ioo 

{t-w)-"dw / s-^e'^^ds 



lOQ 



* a-i,. ^-QJ B{a,l-a) 
1 (a)l (1 — a) 

Our aim is now to explicit, by means of successive steps, the density V2a = 
V2a{u,t) of Ta(t),t > 0: wc first prove that it satisfies a fractional diffusion 
equation of order 2a and, as a consequence, it can be expressed in terms of 
Wright function. Let 

oo f. 



be a Wright function of parameters rj, (3, then we state the following result. 
Theorem 2.3 The fundamental solution to ^ coincides with 

/•oo 

Ua{x,t)^ / Pn{x,u)v2a{u,t)du, (22) 
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where 

V2a{u,t)^^w(^-^;-a,l-ay u > 0, t > 0. (23) 

Proof 

It is proved in Orsingher and Beghin (2004) that, for < a < 1, 

!/■*-- 

r(i - a) Jo 

coincides with the solution V2a{u,t) of the foUowing initial- value problem, for 
< a < 1, 

^v{u, t) = -^v{u, t) u G M, t > 
v{u,0) = 5{u) , .s 

f.(«,0) = ' ^'^^ 

lim|„|^oo v{u,t) = 

where the second initial condition applies only for a G (1/2,1). As a conse- 



quence, formula (|T9)) can be rewritten as ([22)1 with 

- / ,N / 2v2a{u,t), foru>0 , , 

«2«(«,i) = |o, foru<0 • ^2^) 

Since it is known (see, among the others, Mainardi (1996)) that the solution 
to (p4|) can be expressed as 



2t" ^ kir(-ak+l - a) 

= r^T^ I - — ;-a,l- a I , ugM, t>0, 
2i" V J 

we immediately get □ 
Remark 2.1 

By means of the previous result we can remark again that the random time 
Ta possesses a true probability density, which is concentrated on the positive 
half line and moreover it is possible, thanks to representation ([23]) . to evaluate 
the moments of any order (5 > of this distribution. We recall the well known 
expression of the inverse of the Gamma function as integral on the Hankel 
contour 

1 1 



r(a;) 2TTi J Ha 
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which implies the representation of the Wright function as 



k\T [r^k + P) 
— / e > — rfr 

1 / ^-/Jgr+.r-"^^^ 



2ni 



Ha 



Therefore we can show that 

f +C30 



u^V2c,{u,t)du (26) 



/ —W\ ;-a,l-a 



du 



-U^du( ^y-^y^ .-^^ 
t^^^^Jua 



/ eyy-^^-^dy / e-'^-^dz 



r „ . 1 

2TTi 

r(i + (5)r* _ r*r((5) 



r(l + a(5) aT{a5) ' 
From pSjl it is again evident that Jq°° V2a{u,t)du = 1 by choosing (5 = 0. 
Remark 2.2 

It is interesting to analyze the particular case obtained for a — 1/2: in- 
deed, from the previous results, we can show that the process governed by 
gJ/2 u{x, t) = kn-§^u{x,t), a; € R, i > 0, can be represented as {\B{t)\) ,t > 
0, where B{t),t > denotes a standard Brownian motion. This can be seen 

by noting that S'1/2 ('T' -'^''-') coincides with the Levy distribution, so that the 
fractional integral in (I19p reduces to 

r_ I If* we^"^/''"' 

iam[p,M-^;u)\it) = i ^^^^^=^dn. (27) 

u > 0,t > 0, 



TTt 



where the second step follows by applying formula n. 3. 471. 3, p. 384 of Grad- 
shteyn and Rhyzik (1994), for fi — 1/2. Formula ((27)) represents the density 
of a Brownian motion with reflecting barrier in m = 0. This result is con- 
firmed by noting that equation ([24|l . for a = 1/2, reduces to the heat equation 
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§iv{x, t) = -^v{x, t) and then the corresponding process coincides with a Brow- 
nian motion with = 2t. Alternatively, from (j23p . by applying some known 
properties of the Gamma function, we can write 

Vl(u,t) = —=y -^^ TTTT (28) 

_ (-i)fc/2(ut-i/2)fcr(^) 

k even 

_ ^ {-ifl^jut-^'^fT (fc + 1) 0F2i-(^+i) 

^Vt h ^!r(| + i) 



oo 



(-l)JM^J(4t)-J _ e-"'/-^* 

3 On the moments of the solution 

We are now interested in evaluating the moments of the solution of equation 
(HI, that is the moments of the pseudoprocess ^„(7^(i)), i > 0: as we will see, 
they can be obtained in two alternative ways. 

By using the representation of the solution derived in and thanks to the 
independence of the leading process from the (random) temporal argument, we 
can write the r-th order moments as 

E{^l{T^m (29) 



for r g N, t > 0. The moments of the non-fractional n-th order pseudoprocess 
can be evaluated by means of the Fourier transform of the solution of equation 
([2]) which can be expressed as follows 



E 



+00 



(e'^*"W) = / e^^>„(:E,t)da: = e(-''')"'^-"* (30) 

^ / -I — r^ 

(«/3)"^' (-l)"^fc^tJ(nj)! 



— OO 
OO 



Therefore we get 



I r nj 
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which, inserted together with into (P^ . gives, for r — nj, j — 1, 2, ... 

E{KiTa.m (31) 
f-ll"-''/fc-'7n7V f°° 



^ ' " r(aj + 1) 

while it is equal to zero for r ^ nj. 

We can alternatively derive the moments of the pseudoprocesses by evalu- 
ating them directly from the characteristic function of the solution. The latter 
can be obtained by performing successively the Fourier and Laplace transforms 
of equation ([4]) as follows: let us denote by Ua{(3, t) the Fourier transform of the 
solution, i.e. 

r+oo 



then we get form ^ 

^(/3,t) = fc„(-*/3)"^I„(/3,i). (32) 
By applying now the Laplace transform to (15^ we get 
s"C/„(/3,s) - = fc„(-i/3)"C/„(/3,s), 
so that the Fourier-Laplace transform of the solution can be written as 

C/„(^, s) = f / . (33) 

Now recall that for the Mittag-Leffter function 



the Laplace transform (for /5 = 1) is equal to 

e-"'^E„,i(cz")dz 



(see Podlubny (1999), formula (1.80) p. 21, for fc = 0, /3 = 1); hence from 133]) 
we get the following expression for the characteristic function of the solution 

ujp,t)^E^^i{k^{-t(3rn (34) 

and for the solution itself 

u^{x,t)^— e-*-''E„.i(fc„(-z/3)"t")d/3. (35) 
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In the particular case a = 1 the Mittag-Leffler function reduces to the ex- 
ponential so that (|34p coincides with the Fourier transform of the solution to 
the n-th order equation, reported in pop . as it should be in the non-fractional 
case. Analogously, from (P5|) we get the usual expression of reported in 

1^. On the other hand, in the fractional case {a ^ 1) formula fM)) reduces, for 
n = 2, to the well-known Fourier transform of the solution to equation ([3]). 

Finally we can evaluate the moments of the solution by rewriting formula 
dSl as 

(nj)! r(«j + l) 



ii„ (/?,<) - Y. 7r7Tr^i^77-^r(jy + 1), 



so that we get again expression ([3T 



4 More explicit forms of the solution 

In order to obtain a more explicit form of the solution to , in terms of known 
densities, we need to distinguish between two intervals of values for a. 

(i) Case 1/2 < a < 1 

If we restrict ourselves to the case a G [1/2, 1) , so that 1 < 2a < 2, it is 
possible to apply a result obtained in Fujita (1990), which expresses the solution 
to a time- fractional diffusion equation in terms of a stable density of appropriate 
index. By adapting that result to our case, we can conclude that the solution 
to ([M)) coincides with 

where Pi/q,(-; i) denotes a stable density of index 1/a G [1,2) with parameters 
a = (fcos(7r - ^))", /3 = -1, = (for brevity 5i/„((7, -1,0)). 

Therefore the density of 7^(t), i > is proportional to the positive branch 
of a stable density, as the following expression shows: 

W2q(w, i) = — Pi/q(u; i), u > 0, i > 0. (36) 



Remark 4.1 

It is possible to recognize, in the previous expression, a known density, by 
resorting to results on the supremum of stable processes (see, for example, Bing- 
ham (1973)). More precisely, let us define Y{t) — supg<s<( Xi/(j(s) where 
Xijait), i > is a stable process of index 1/ a and with characteristic function 

^(e-^i/.(t)) =exp|-t|s|i/" (^1-i-itan^^^l , t,s>Q. 

It corresponds, for any fixed i, to the stable law Pi/a{'',t) defined above 
and, for t varying, to a spectrally negative process, which has no positive jumps 
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(since, for /3 = —1, the Levy-Khinchine measure assigns zero mass to (0,oo), 
see Samorodnitsky and Taqqu (1994), p. 6). Under these circumstances and for 
1/a G [1,2) , it is known that the Laplace transform of Y{t) is equal, for any 

s,t > 0, to 

i?(e-^^(*))=E„,i(-si"), 

where Yia,i3{x) is the Mittag-LefHcr function defined above. Since it is also well- 
known that 



e ''^pi/aiu;t)du ^ aEa,i{-st"), t,s>0, 
we can conclude that 



a 



Alternatively it can be shown, by adapting the result of Bingham (1973), 
that the density of Y{t) can be written as 

P{Y(t)edu} = > ^- sin(7rna)r(l + na) — du 

an ^-^ n\ Vi"/ 

n— 1 

— —Pi/ai'U',t)du, u > 0, t > 
a ' 



which coincides with ([36]) . 

Formula ([55)1 shows that, for 1/2 < a < 1, 

lil-a) [Pa{w;u)] (t) ^ -Pi/a{u;t), U>0,t> 0. 

a 

Then, as a result of the fractional integration of the stable density p^{-;t), 
which is totally skewed to the right (with support [0, oo)), we obtain the positive 
(normalized) branch of a new stable density Pi/ai-'it) (defined on the whole real 
axes, since it is 1/a € (1, 2]), which represents the distribution of the maximum 
of a stable process of index 1/a. 

(ii) Case < a < 1/2 

We turn now to the other interval of values for a, i.e. (0, 1 /2] , so that, in this 
case, it is < 2a < 1. An explicit expression of the solution can be evaluated 
by specifying a = 1/m, m S N, m > 2. In this particular setting, problem (I24|) 
becomes 

^v{u, t) = -^v{u, t), u e M, t > 
v{u,Q)^5{u) (37) 
lim|„|_oo v{u,t) = 

so that it can be considered as a special case of the fractional telegraph equation 
studied in Beghin and Orsingher (2003), for A = and c — 1. By applying 
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formula (2.11) of that paper, the solution to ([57| can be expressed, for u G M, 
t > 0, as 



V2/ni{u,t) (38) 

m \ 1 
2^^ io 



oo 

dwi... I dw,yi-i 



• e ^^~^StS^W2 • • • [(5(U - Wi • • • Wm-l) + S{u + Wi ■ ■ ■ Wm-l] ■ 

By taking, as before, 

w,/ - / ^''2/™(w,i), forii>0 , . 

V2/m[U,t)-<^ 0, forii<0 ' ^"^^^ 

the solution to our problem ^ can be expressed, in this case, as 

/>oo 

Ul/m{x,t) = j Pn{x,u)v2/m{u,t)du 

) 

Pn{x,u)pG(t){u)du 



where G{t) — Ojli^ G'j(i),i > and Gj{t), j = 1, ...,r7i — 1 are independent 
random variables with the following probability law 

1 / \ 

PG,(t){w) = , exp - rr^-ir-^ ^'^^ w>0. (40) 

We can check the independence by noting that 

m — 1 

n PGAt)iw3) (41) 



m — 1 ^ 
1 J. 7 1 J exp Tj, — 1 y 

exp ( - 1 , 



27r/ ^ \ "^-y/m^t 



where, in the second step, we have applied the multiplication formula of the 
Gamma function 

r(z)r (z + —\ ...r (z + ^ (2^)^ m5-"-r {mz) , 

\ m J \ m J 

for z = 1/m. The last expression in (|4ip coincides with the joint density of the 
variables Gj{t) given in formula (1.7) of Beghin and Orsingher (2003). Therefore 
the corresponding pseudoprocess is represented, in this case, as "^n (G(<)) , t > 0. 
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Remark 4.2 

We can check the previous results, obtained separately for the two intervals, 
by choosing a = 1/2. From both cases we obtain again that the pseudoprocess 
governed by our equation can be represented by , t > 0. 

Indeed from the first case, i.e. for 1/2 < a < 1, we get, by means of 
that the density oiTa{t),t > 0, for a — 1/2, coincides with the folded normal. 
More precisely, S'2(\/i, — 1, 0) coincides with N{0,2t) and then 

v,(u,t)^2p2{u;t) = ^^ (42) 
Vnt 

for M > 0, t > 0. 

On the other hand, if we consider the expression of the density oiTa obtained 
for < a < 1/2, we get, for a — 1/2 and m — 2, from (|¥T|) that again 

PGi(t)(w) = — (43) 
Moreover both ([1^ and coincide with ^7} derived above, as expected. 

An interesting application of our results can be obtained by specializing The- 
orems 2.2 and 2.3 to the particular case n = 2. In this situation the pseudopro- 
cess 4'„(i),t > reduces to the Brownian motion (with variance 2t) B{t),t > 
and therefore the solution of Q coincides with the transition density of the 
process B{Ta{t)), t > obtained by the composition of B with the random time 
Ta- We state this last result as follows 



Corollary 4.1 

The solution of the problem 

^u(x, t) = £iu{x, t) a; e M, t > 0, 
u{x,f))=5{x) ^ ' 

for < a < 1, is represented by the transition function of B{Ta). The density of 
the random time Ta ^ Ta{t) ,t > Q \s the folded solution of the time- fractional 
equation 

and is given in (j23|) . 



We can prove that this is in accordance with what is already known on (|44|) : 
for n = 2 we can substitute in (|22|) the transition function of the Brownian 
motion, so that we get: 
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If we prove now that 



ey , e~^-^y , 

ay / du. 







-du = y^r/2y-"/2e-I^^ (46) 



and substitute (HS)) into we finally get the known result 



1 1 f 



1 W/'. N « « 



2t"/2 V t"/^' 2' 2 



In order to verify formula (j46[) we use the following relationship, known for 
the Laplace transform of the first-passage time of the Brownian motion, 

Jo 2y^V 

which, integrated with respect to x gives (|46|) . for s = Alternatively, we 

can apply formula n. 3. 471. 9, p. 384 of Gradshteyn and Ryzhik (1994), for /3 = 
x^/4, 7 = = 1/2 (noting that ii'i/2(-z) = ^7r/2ze~^, see Gradshteyn 

and Ryzhik (1994), n.8469.3, p.978). 

If we restrict ourselves to the case a G [1/2, 1) , the density of Ta{t),t > 
is again proportional to the positive branch of a stable density, as expressed in 

On the other hand, for a e (0, 1/2] and in particular a — \/m, m € N, it is 
represented by the law presented in and ([55]) , so that the process governed 
by equation (|44|) is, in this case, i^dlj'Ti^ Gj(^)): ^ > 0. 
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